Elliptic gamma-function and multi-spin solutions of the 

Yang- Baxter equation 



Vladimir V. Bazhanov 

Department of Theoretical Physics, 
Research School of Physics and Engineering, 
Australian National University, Canberra, ACT 0200, Australia. 

and 

Mathematical Sciences Institute, 
Australian National University, Canberra, ACT 0200, Australia 
E-mail: Vladimir.Bazhanov@anu.edu.au 

Sergey M. Sergeev 

Faculty of Information Sciences and Engineering, 
University of Canberra, Bruce ACT 2601, Australia. 
E-mail: Sergey. Sergeev@canberra. edu. au 



Abstract 

We present a generalization of the master solution to the quantum Yang-Baxter equation 
(obtained recently in arXiv:1006.0651) to the case of multi-component continuous spin variables 
taking values on a circle. The Boltzmann weights are expressed in terms of the elliptic gamma- 
function. The associated solvable lattice model admits various equivalent descriptions, including 
an interaction-round-a-face formulation with positive Boltzmann weights. In the quasi-classical 
limit the model leads to a new series of classical discrete integrable equations on planar graphs. 
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1 Introduction 



There are solvable lattice models of statistical mechanics with only a pair interaction between 
neighbouring spins, i.e, where two spins interact only if they are connected by an edge of the 
lattice. The Yang-Baxter equation for these models usually takes the form of the "star-triangle 
relation" [lj. The most notable discrete-spin models in this class include the Kashiwara-Miwa [2] and 
chiral Potts [3H2] models (both of them also contain the Ising model pQ and Fateev-Zamolodchikov 
Z/v-model [6] as particular cases) see [7j for a review. There are also important continuous spin 
models, including Zamolodchikov's "fishing-net" model [8j, which describes certain planar Feynman 
diagrams in quantum field theory, and the Faddeev-Volkov model [9], connected with quantization 
[10] of discrete conformal transformations [111112] . 

Recently [13] we have found a new solution of the star-triangle relation which contains as special 
cases all the discrete- and continuous-spin solutions mentioned above3, and also leads to new ones. 
This "master solution" is expressed through elliptic gamma-functions and contains two temperature- 
like variables. It defines an exactly solvable lattice models with continuous spin variables taking 
values on a circle. Its connection to the theory of elliptic hyper geometric functions is discussed 
in |14j . From an algebraic point of view the model is related to the modular double [151116] of 
the Sklyanin algebra [17] . The latter is an elliptic deformation of the quantum group U q (sl(2)), 
connected with the i?-matrix of the eight- vertex model [18] . 

In this paper we extend the main results of [13J to the case related to the sl(n) algebra with 
n > 3. The generalized model has multi-component continuous spin variables taking values on a 
circle. Similarly to the n = 2 case of [13] the model contains two temperature-like parameters. The 
Boltzmann weights satisfy the so-called star-star relation pL9l|2Q] (see Eq. (|26p below), which ensures 
the integrability of the model. Currently, we claim this relation as a conjecture, however, we expect 
that a complete proof could be obtained by a rather straightforward generalization of the results 
of [201424] devoted to discrete-spin models connected with the sl(n) algebra. 

It should be noted, that apart from the n = 2 case previously considered in [13], the two-spin 
Boltzmann weights are not real and positive. Fortunately, this is not an indication that the model 
is unphysical. It can be reformulated as an "interaction-round-a-face" (IRF) model and then there 
exist a domain of parameters, where the IRF-type Boltzmann weights become real and positive, see 
Section [3j 

The quasi-classical (or low-temperature) limit of the model is considered in Section [H A sta- 
tionary spin configuration which gives the leading contribution to the partition function in this 
limit is described by new classical discrete integrable equations for multi-component fields assigned 
to lattice sites. These equations can be thought as a generalization of the Laplace-type equation 
associated with the famous system [25, 26j to the multi-component case. 

In Conclusion we summarize the results and discuss their connections with some other integrable 
systems in two and three dimensions. 

2 Formulation of the model 

We start with the definition of the elliptic gamma- function |16p27ff29] . Let q , p be two elliptic nomes 
(they play the role of the temperature-like parameters), 

p = e ,na , q = e i7rT , 1m a > 0, lmr>0. (1) 

To be more precise, it only contains the solutions, which have a single one-dimensional spin at each lattice site. 
For this reason, it cannot contain the D > 2 fishing-net model which has multi-dimensional spins. 
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In principle, these parameters can be arbitrary (apart from the requirements |p| < 1 and |q| < 1), 
however, in the following we will often refer to special regimes when p and q are either real or 
complex conjugate to each other, 

(i) p* = p, q*=q, (ii) p*=q. (2) 

Note that in both of these cases the "crossing parameter" 

77 = -w(ff + r)/2, (3) 

is real and positive. Define the elliptic gamma-functiorj! 

00 1 - e 2i ' z q 2 J+ 1 p 2A ''+ 1 ( e -2\zk | 

= II i_ e -2i 2q 2i+i p 2fc+i = ex PS Efc(q*- q -*)(pt- p -*) f ' ^ 

where the product formula is valid for all z, while the exponential formula is only valid in the strip 

— Re r\ < Im z < Re rj . (5) 
The function @ possesses simple periodicity and "reflection" properties 

$(z + vr) =$0), $(*)$(-*) = 1. (6) 
Moreover, it satisfies the following difference equation 

l^'tl = n( 1 - e2i V n+1 )( 1 - e_2i V"' +1 )=^|r), (7) 

^\ Z+ 2 ) „=0 

and a similar equation obtained by interchanging r and a. Here 

^kH^-y^lr), G(r) = [I(l-q a ), q = e iwT , (8) 
^ ' fc=i 

where #j(2|t), j = 1,2,3,4, stand for the standard elliptic theta-functions of the periods tt and 
7tt, as defined in [30]. Note also, than in the regimes ([2]) the function @ has a simple complex 
conjugation property, 

$(*)* = *(-**). (9) 



Next, we want to introduce a new two-dimensional solvable edge-interaction model. The model 
can be formulated on rather general planar graphs, however, for the purposes of this presentation 
it is convenient to take a regular square lattice. Consider the square lattice, drawn diagonally as in 
Fig. [TJ The edges of the lattice are shown with bold lines and the sites are shown with either open 
or filled circles in a checkerboard order. In this Section we will not distinguish the two type sites; 
their difference will be important in Sect. El At each lattice site place a n-component continuous 
spin variable 

n 

x = {xi, . . . ,x n } G R n , 0<Xj<TT, ^ Xj = (mod 7r) . (10) 

i=i 

2 Our function &(z) coincides with r(e~ 2l ' z_ '' , p 2 , q 2 ) in the notation of ref. |16| . The definition Q differs slightly 
from that of ref. [T3], where the RHS of @ is denoted as <E>(2z). 
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Figure 1: The square lattice shown with bold sites and bold edges drawn diagonally. The associated medial 
lattice is drawn with thin and dotted horizontal and vertical lines. The lines are oriented and carry rapidity 

variables u, u' ', v and v'. 



Note that due to the restriction on the total sum, there are only (n 
For further reference define the integration measure 



1) independent variables Xj. 



dx 



dx\ ■ ■ ■ dx n -i , 1 



\x)dx{x\ 



(11) 



Fig. [T] also shows an auxiliary medial graph whose sites lie on the edges of the original square 
lattice. The medial graph is drawn with alternating thin and dotted lines. The lines are directed 
as indicated by arrows. To each horizontal (vertical) line on the medial graph assign a rapidity 
variable u (v). In general these variables may be different for different lines. However, a convenient 
level of generality that we shall use here is to assign the same rapidity u to all thin horizontal lines 
and the same variable v! to all dotted horizontal lines. Similarly, assign the variables v and v' to 
thin and dotted vertical lines as indicated in Fig. [H 






U X V 



- v >(y, x ) 



Figure 2: Four different types of edges and their Boltzmann weights. 



Two spins interact only if they are connected with an edge. The corresponding Boltzmann weight 
depends on spins at the ends of the edge and on two rapidities passing through the edge. There are 
four types of edges differing by orientations and types of the directed rapidity lines passing through 
the edge. They are assigned with different Boltzmann weights as shown in Figj2j The weights are 
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defined as 

n 

W a (x,y) = K n (a)- 1 H $( Xj -y k + \a), W a (x,y) = ^/§(x)§(y)W v - a (x,y) , (12) 
j,k=i 

where the single-spin function § is given by 

S(sc) = ^^{^xj - Xk + ^y 1 = ^ll^Mx, - Xklr^^xj - x k \a)^ , (13) 

j^k j<k 

and 

- - - ! (o^r • <i4) 

Here the indices j, k run over the values 1,2, ... ,n; the functions $(x | r) and G(t) are defined in 
([8]). The normalization factor 

^ (p fe - p" fe )(q fe - q _fe ) p nfc q nfc - p -nfcq-"fc j ' ( 15 ) 

has the meaning of the partition function per edge for unnormalized Boltzmann weights (i.e., when 
the factor K n (a) in (|12[) is omitted). It solves a pair of functional equations 

K n {a) n n (—a) = 1, K n (r] — a) n n (r] + a) = <3?(ir7 — \na) ^(ir/ + ma) . (16) 

The weights (|12p satisfy two inversion relations: 

W a (x,y)W- a (y,x) = 1 , J dxW a (x, z)W^ a (z,y) = ^Y, 6 ( x ^^ ( 17 ) 

where the sum is taken over n! permutations <r of components of the vector y = {yi, . . . ,y n } and 
the 5- function is understood with respect to the measure (jlip . The first of these relations is a trivial 
corollary of the definition f 1 1 2 [) and the reflection property of the elliptic gamma function &(z). The 
second relation is a particular case of a more general relation stated in Theorem 11 in ref. |16j . 
Note that for any permutation a one has 

W(aj, y) = W a (a(x),y) = W(x, &(y)) , S(x) = S{a(x)) , (18) 



as a trivial consequence of the definitions (I12p and f)13|) . Under the complex conjugation the weights 
transform as 

W a (x,y)* =W a *(y,x)=W a *(-x,-y) , S(x)* = S{x). (19) 

provided the nomes p and q belong to either of the regimes, defined in ([2]) (remind that the spin 
variables are always assumed to be real; see (|10[> ). When n = 2 the weights W and § are real 
and positive [T5] . Correspondingly, the weights W a (x,y) and W a (x,y) are symmetric under an 
exchange of the spins x and y. However, for n > 2 this symmetry is lost. As a particular consequence 
of this fact the above description of the associated lattice model required two types of the rapidity 
lines. 

The square lattice in Fig. [T]can be formed by periodic translations of the "box diagram" shown 
in Fig. [3l The Boltzmann weight of this box can be conveniently associated with an i?-matrix, 

(x,y\R uv \x\y') =W u - v (y,x')W u '- v/ (y',x)W u ^ v (x',y')W u ^(x,y) , (20) 



5 




where u = [u, u'\ and q = [v, v'] stand for the pairs of rapidities in horizontal and vertical directions. 
The partition function is defined as 

Z = j [] (x,y\^uv\x',y') l[dx, (21) 

boxes sites 

where the first product is taken over all boxes and x,y,x',y' are the spins at the corners arranged 
as in Fig. [3j The integral is taken over all configurations of the spin variables on the internal 
lattice sites. The boundary spins are kept fixed. Note that due to periodicity ([6]) the definition 
(I2ip only contains closed contour integrals (the contours can be deformed into the complex plane, 
if necessary). 

3 Yang-Baxter equations and star-star relations 

The i?-matrix (|20p can be regarded as a kernel of an integral operator acting on a pair of (multi- 
component) continuous spin variables (|10p . We claim that it satisfies the Yang-Baxter equation of 
the form 

j dx'dy'dz'(x, y \R uv \x', y') (x', z \R uw \x", z') (y' , z'\R vw \y", z") = 

(22) 

J dx'dy'dz'{y, z\R vw \y', z') (x, z' \R uw \x' , z") (x' , y'\R uv \x", y") , 

where the integration measure is defined in (jlip . and the symbols u = [u,u'], v = [v,v'\ and 
w = [w, w'] stand for the rapidity pairs. With the standard conventions the last equation can be 
written in an operator form 

KLuu M vw = IRitio Ruij • (^3) 

The i?-matrix (|20p was derived as a unique intertwiner for two different sets of Lax operators serving 
the modular double [T5l[T6] of the sl(n) analog of the quadratic Sklyanin algebra [TT|l3T]. The details 
of calculations will be published elsewhere. 

According to the terminology of the Baxter's book |32j . Eq. (l2ip defines a vertex model. It is 
easy to see that to within boundary effects the same model can be equivalently reformulated as an 
interaction-round- a-f ace (IRF) model. Indeed, the lattice in Fig. [T]can also be formed by periodic 
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translations of a four-edge star (consisting of four edges meeting at the same site), instead of the 
box diagram of Fig. [3l A little inspection shows that there are only two different types of such 
stars shown in Fig. HI Recall that there are types of sites, shown with open and filled circles. We 
will refer to them as to "white" or "black" sites, respectively. There are white-centred stars (i.e., 
centred around a white sites) and black-centred ones. Applying the rules shown in Fig. [2] one can 
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Figure 4: Two types of four-edge stars: a white-centred star V^ 1 ) (left) and a black-centred V^ 2 ^ (right) 



write IRF-type Boltzmann weights corresponding to these stars 

dx W u - V (c, x) W u /-„' (6, x) W u >- V (x, a) W u - V > (a;, d) , 



(i) / a b 
Vuv[ c d 



and 



V 



(2) 



a b 
c d 



dy (y,b) (d,y) 



(24) 



(25) 



where the bold symbols u and v has the same meaning as in (1221) . The above two expressions are 
connected by the so-called star-star relation [19, 20J , which in our case reads 



r v ,_ v (d,c)w u ^ u (d,b) v[ 



(i) f a b 
c d 



7 V '- V (b, a) W u /_ u (c,a) Y^l ( ° ^ 



(26) 



Apparently this is the simplest relation for the Boltzmann weights which ensures the integrability 
of the considered mode@. Currently, we claim this relation as a conjecture, however, we expect 
that a complete proof could be obtained by a rather straightforward generalization of the results 
of [20H24"] . devoted to discrete-spin models connected with the sl(n) algebra. We have verified this 
relation in a few orders of perturbation theory in the temperature-like variables (see Appendix |B|) 
and made extensive numerical checks for n = 3,4. 

It is worth noting that the Yang-Baxter equation (|22p for the composite "box i?-matrix" (|20p 
is a simple corollary of the star-star relation (I26p . The same relation also implies yet another 
Yang-Baxter equation for the IRF-type weights 



dh V.. 



c h 
e d 



h b 
d f 



da V, 



"trai 

e a 
d f 



c g 
h b 

Vim 



c g 

e a 



V,. 



g b 

a f 



(27) 



3 For n = 2 the star-star relation (|26p is just a consequence of the star-triangle relation, Eq.(1.5) of [T5], which is 
equivalent to the elliptic beta integral [141129]. However, for n > 3 the corresponding star-triangle relation apparently 
does not exist (at least it is not known to the authors) and the star-star relation (|26[) seems to be the simplest relation 
of this type. 
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where 



, a b \ J W v >- v (d,c)W u >- u (d, b) \ 2 v (i) fab 

c d J " 1 W„,_„(&,a)W u ,_ u (c,a) f c d 



i 



V_„(b,a)W u /_ u (c,a) 1 2 v (2) / a 6 
J„,_„(d,c)W u ,_ u (d ) &) / Vu M c d 



i 



(28) 



The partition function (12ip can be re-written using either of the weights (|24l) . (j25|) or (128j) . For 
example, it easy to see that, up to boundary effects, 

z= j n v «« ( " d ) n ( 29 ) 

(white stars) " (black sites) 

where the first product is taken over all white-centred stars and a, b, c, d denote the corresponding 
outer spins, arranged as in Fig. 01 The integral is taken over the spin variables on black internal 
sites (the integration over the spins on white sites is included through the definitions (|28p and (|24D ) . 
As before, the boundary spins are kept fixed. Note that the weights (124H . (125p and (128|) differ from 
each other by equivalence transformations which leave the partition function (j29[) unchanged (up 
to boundary contributions). 

There are two alternative forms for the RHS in (|28p . which coincide by virtue of the star-star 
relation (j26|) . This fact allows one to understand the behaviour of the IRF- weight (j28|) under the 
complex conjugation. Consider any of the regimes defined in ([2]) and assume that the rapidity 
variables obey the relations 

u* = u , v* = v' , and < Re(u' — v' ,u — v, . . .) < rj . (30) 

Note that these relations allow a homogeneous case u = u', v = v', when u and v are real. Taking 
(f30|) into account and using (fT9|) one can easily check that the complex conjugation just interchanges 
the two equivalent forms of (128|) . It follows then that the IRF weight (128p is real. Further, one can 
show that there must exist a region for the parameters p and q where this weight is non-negative, 



a b 
c d 



V m ,( I b d ]>(.). (31.) 



for all values of the spin variables. First of all note, that in the regimes ([2]) the single-spin weight 
(|13p is real and non-negative 

S(£c)>0. (32) 

It vanishes only then at least two components of x coincide. Next, the definition (|12p for W Q (cc, y) 
contains square roots of S(x) and S(y) and therefore there is an ambiguity in choosing signs. 
However, it is easy to see that these signs cancel out for all internal sites and therefore can be 
chosen arbitrarily. For definiteness we assume that S(:c)2 = |S(sc) 2 | > 0. 

Consider the limit when p, q — > and the ratio p/q is finite. It is convenient to parametrize the 
rapidity variables as 

« = 5+0 + 5(7-/8), v = -~(7 + /3), 

: (33) 

«' = \ +«" ^(7-/3), v' = +^( 7 + /3) 
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where the variables a, (3 and 7 are real. Moreover, assume that 

|Rea|< 77, (34) 

which means that e 2a = 0(1), when p,q — > 0. The relations (|30p are automatically satisfied. The 
weights (fT2j) and (fT3|) can be expanded in powers of p and q 

W a (x,y) = l + 0( P q) + ..., 



W v/2+a (x,y) = l + 0((pq)2) + O(pq) + 

n^nl S(x) Yl (2 sm( Xj - x k ))~ 2 = 1 + 0(p 2 + q 2 ) + . . . . 

j<k 



(35) 



Explicit form of the coefficients to within the forth order is given in Appendix [B] (the fractional 
powers arise due to the relation e~ v = (pq)a ). It is convenient to define sums of exponents of the 
spin variables, entering the IRF weight (|28p 



A 1 = Y^^, £i=^e 2i \ Ci = ^e 2ic ^, X>i = ^e 2i ^, (36) 
j=i j=\ j=i j=i 



Using the above expansions one obtains 



S(6)-fS(c)-f Y uv [ a b . ) =l + pqP + 0((pq)t) + ... (37) 



c d 

where V > 0, denotes the following expression 

V = (e 2i/3 Ax + e" 2i/3 Pi) (e 2i7 B{ + e" 2 ' 7 C{) + (e" 2i/3 A\ + e 2i/3 V{) (e" 2 ' 7 B x + e 2 ' 7 Ci) , (38) 

which depends on spins and rapidity variables. A few more terms of the expansion is given in 
Appendix |Bj They are all manifestly real and thereby confirm the validity of the star-star relation 
(I26p in perturbation theory. Note that the leading term in (I37p is strictly positive. The coefficients 
of the expansions in (|35p are analytic functions of the spin variables, thus the integrals over central 
spins in the star weight (I28p are non-singular in every order of the expansion. It follows then that 
the coefficients in front of powers of p and q in the RHS of (|37|) are analytic function of the external 
spins. Therefore, there exists a finite domain of the parameters p and q in the vicinity of p = q = 0, 
where the RHS of (|37p is strictly positive for all values of the external spins a, b, c, d. We will call 
this domain a physical regime. It would certainly be interesting to investigate exact boundaries of 
this domain. 

In the physical regime the considered lattice model becomes a well defined model of statistical 
mechanics and Euclidean lattice field theory. Using the standard arguments based on the commu- 
tativity of transfer matrices and functional equations for the partition function per edge [33H35] . 
we deduced that for a large number of sites N (in the thermodynamic limit) the bulk free energy 
of the model vanishes, 

lim N^ 1 log Z = , (39) 

N^oo 

provided that the normalization (I15p in (|12p is taken into account. 
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4 Quasi-classical limit 

4.1 Asymptotics of the Boltzmann weights 

Let us now consider the quasi-classical limit of our model, when one of the nomes is real and fixed, 
while the other one tends to unity, 

q = e ' nT = real, p = e iwa 1, a -)• . (40) 

It convenient to use a rescaled parameter 

n=-\Tra^0, (41) 

which plays the role of the Planck constant from the point of view of Euclidean field theory (or the 
temperature from the point of view of classical statistical mechanics). Introduce a new function 

\ a (z\t) = -i f dw logd A (z \ t) , (42) 
Jo 

where I &j(z\r) is defined in ([8]). In the limit (140p the elliptic gamma- function (J4]) and the factor 
(fT5l) become 

log$(*) = --M{z\T)/h + O(h ), log Kn (a) = -^X 4 (ma\nT) + O(h ), (43) 

moreover 

logG(a) = -^-hogh+Oih ), \ogMx\a) = ~{x 2 -7rx+^}+O(h ), < x < vr , (44) 

where G{a) is defined in ©. In writing quasi-classical expansions of the Boltzmann weights it is 
convenient to employ the permutation symmetry (1181) . Below we will assume that components of 
spin variables are always arranged such that 

7T 7T 

x = (x 1 ,x 2 ,...,x n ), -— < Re(xi) < Re(x 2 ) < • • • < Re(x n ) < - . (45) 
The leading quasi-classical asymptotics of the weights (|12p reads 

\ogW a (x,y) = -lA a (x,y) + O(h ), (46) 
n 

logW a (x,y) = -^A a (x,y)-^^\ogh+O(h ), (47) 



where 



A a (x,y) = —\^(\na \ nr) + \i{xj — yp, + la \ r) , (48) 



A a (x,y) = A V0 - a (x,y) + ^C(x) + ^C(y), rjo = -ivrr/2. (49) 

Eqs.(|48p and (|49p are trivial corollaries of the definition (|12p and the expansions ()43p . The function 
C(a;) is determined by the asymptotics of the single-spin weight 

log §(*) = ~c{x) - ^^\ogh +o(h°), (50) 
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defined in (fT3"j) . Using (|44p one obtains 

( n 2_ l \ n 2 n n 

C(x) = ± ^— + n ~ 2?r J *i • (51) 

3=1 3=1 

provided the variables Xj are arranged as in (|45p . Finally, taking into account that q is assumed 
real, one deduces that under the complex conjugation 



A a (x,y)) =A a *(y*,x*), (52) 
where x* = (a^, x|, . . . , x*) and similarly for y* . 

4.2 Classical star-star relation 

Consider now the quasi-classical limit for the composite Boltzmann weights of the stars shown in 
Fig. [U To make equations more compact we will write the spin arguments a, b, c, d on the same line, 
but assume the same spacial arrangement as in Fig. |H Substituting (|4"T|) into ([24"]) and calculating 
the integral by the saddle point method, one gets 

log (V£2(a, b, c, dj) = -~ 6uU X I a > b > c > d ) ~ log H + > (53) 



where 

(as | a,b,c,d) = A u _ v (c,x) + A u /_ v /(6,£c) + A u /_„(a;, a) + A u _v(a;,d) (54) 
and X = (Xi,X2, ■ ■ ■ ,X n ) is the saddle point, determined by the equations 

(n * ) C$ v (x\a,b,c,d) =0, fc = l,2,...,n-l. (55) 

\OXk OXfc+l/ x=X 

Remind that x is the multi-component variable, x = (x% , x-i , ■ ■ ■ , x n ) , whose components x k are 
constrained by (|1U|) and cannot be varied independently. That is why there are only (n — 1) 
equations in f|55|) . each containing a difference of two partial derivatives. Define new variables 



(56) 



a% = v! — v , «2 = 770 — v! + v ' , as = 770 — u + v , 

«4 = U — V 1 , 05 = u' — U , «6 = w' — V , 

constrained by three relation 

a\ + a 2 + 03 + 04 = 2770 , 05 = «i + 03 - ry , «6 = + «2 — % • (57) 
Below we will regard 01,02,03 as independent variables. Introduce a new function 



1 1 1 , ■ V^X - a k + iai |r)# 4 (X - 4 + ia 4 t) 

# 4 (c fc - ^ + ia 3 |r)^ 4 (ofc - ^ + i«2 |r) 

which depend on a (single-component) variable X and four multi-component spin variables a, b, c, d 
and the rapidity- type variables ai,Q2,«3,a 4 (it is also implicitly depends on the elliptic modular 
parameter r). Using (|48f) . (|49p and (|5ip one can write the variational equations (|55p in the form 

I o, 6, c, d I ai, «2) 03, a 4 ) — ip(X k \ a,b,c,d \ a\, 02, 03, a 4 ) = 

(59) 

= +2vr - 2n (X fc+1 - X k ) , k = 1, 2, . . . , n - 1 , 
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where the RHS comes from a variation of the function C{x) associated with the central spin in the 
star (it enters the expression (|49p ). Similarly, for the second star in Fig. U one gets 

log (v$(a, b, c, d))=~ 6Sv{Y | a, b, c, d) - log h + 0(H°) , (60) 

where 

£u}v (y I a, b, c, d) = A u _„(y,&) + A u ,_ v >{y,c) + A u <_ v (d,y) + A u _„/(a,y), (61) 
and the saddle point V = (Yi, Y2, ... , Y n ) is determined by the variational equations^] 

V'CVfc+i \ a,b,c,d \ -0*4, -a 3 , -a 2 , -«i) - ip(Y k \ a,b,c,d \ -a 4 , -a 3 , -a 2 , -a{) = 

(62) 

= -2vr + 2n (y fe+1 - Y k ) , k = 1, 2, . . . , n - 1 . 

The quasi-classical expansion of the quantum star-star relation (j26j) leads an infinite number of 
non-trivial relations — one relation for each order of the expansion in h. In the leading order one 
obtains, 

Cu]v (X\a, b, c, d) - A u>v (a, b, c, d) = Cujv (Y\a, b, c, d) + A U)V (a, b, c, d) , (63) 
where X and Y are determined by (|59[) and (|62p and 

A u<v (a, b, c,d) = - (A„/_ v (6, a) + A u /_ u (c, a) - A v ^ v (d, c) - A u /_ u (d, b)) . (64) 

We name the formula (163p a classical star-star relation. 
For the weight function (|28|) one obtains 

1 (n — 1) 

logV u „(a,b,c,d) = — -£ u l> (a,b,c,d) — log ft + £>(£) , (65) 

n 2, 

where the Lagrangian density C UtV (a,b,c,d) coincides with the LHS (or the RHS) of the (|63|) . 
Writing it in full, one obtains 

C UjV (a,b,c,d) = A u ^ v (c,X) + A u ,_ v ,(b,X) + A u i_ v (X,a) + A u _ v , (X , d) - A u , v (a,b,c,d) 
= A u _ (y, b) + A U '_ V '(Y", c) + A u /_ v (d,Y) + A 

u—v 

(a, b, c, d) , 
(66) 

where X and Y" are determined by the variational equations (|59p and (|62p . 

Consider now the regime ()3Q|) with real q as in (I40p . Assume that the spins a,b,c,d are real. 
Then one can show that if X solves (|59p then Y = X* solves (|62p . where the star denotes the 
complex conjugation. Having this in mind and using (152p it is easy to check that the complex 
conjugation just interchanges two alternative expressions in (|66p . This means that the Lagrangian 
density C u v (a,b,c,d) is real. 

From ((I2D and flUD it follows that 

A a {a,b) + A_ Q (6,a) = 0. (67) 
4 Notice different signs in the RHS of this equation with respect to (1591) . 
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h c 

Figure 5: Graphical representation of the classical star-star relation (|68|) . 



Using this equality together with (j49j) . substituting (|54|) . (f6Tj) and (j64l) into the classical star-star 
relation (|63p and moving all terms there to one side, one obtains 



Aa! (X, a) + A Q2 (5, X) + A Q3 (c, X) + A a4 (X, d) + A a6 (d, c) + A a5 (d, 6) 
+ A_ a4 (Y , a) + A_ a3 (b, Y) + A_ Q2 (c, y ) + A_ Ql (y , d) + A_ Q6 (a, 6) + A_ Q5 (a, c) (68) 
+ C(X) - C(Y) = 0. 



The terms of this relation can be conveniently associated with edges and vertices of an octahedron. 
Consider a perfect octahedron shown in Fig. [5j Its edges are oriented, as indicated by arrows, 
such that there are exactly two incoming and two outgoing edges at each vertex. The vertices are 
assigned to the spin variables a, b, c, d, X , Y, as indicated. The edges carry rapidity-type variables 
±cki, ±ci2, • • • , i«6> shown near edges. A perfect octahedron contains six pairs of parallel edges. 
The values of the a-variables on parallel edges are equal in absolute value, but differ in signs. 
Moreover, the a-variables are constrained by the relations (157p . Each edge corresponds to a A-term 
in ()68|) . For example, the edge directed from X to a and carrying the rapidity label a 4 in Fig. [5] 
corresponds to the term A Qil (X, a) in ()68p . Similarly for all other edges. Further, according to 
(I57p . the arithmetic sum of the a-variables on four edges meeting at the vertex X (regardless their 
directions) is equal to +2tt. This vertex contributes the term "+C(X)" to ([68^ . Similarly, the sum 
of a's at the vertex Y is equal to —2ir. This vertex contributes the term "— C(Y)" . Finally, thanks 
to (|57p the sum of a's for each of the remaining four vertices a, 6, c, d exactly vanishes. These 
vertices do not contribute any C-terms into (|68p . 

From mathematical point of view it is convenient to regard equations (|59p and (|62p as constrains 
placed on six spin variables at the vertices of the octahedron. Interestingly, these constrains can be 
re-written in two other equivalent forms. Differentiating (j68|) with respect to the spins a and d one 
obtains 



rp(a k+ i | b, X, y, c | -a 6 , +ai, -a 4 , -a 5 ) 



ip(a k | b, X, Y , c | -a 6 , +ai, -a 4 , -a 5 ) = , 



(69) 



ip(dk+i \ b,X,Y,c\ +a 5 , +a 4 , -ai, +a 6 ) 



V>(<4 I 6, X, y , c I +a 5 , +a 4 , — ai, +a 6 ) = 
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where k = 1, 2, . . . , n — 1. Here we have used the fact that the LHS of (|68|) is stationary with respect 
to variations of X and Y by virtue of ([59]) and ([62]) . Similarly, differentiating ([68]) with respect to 
the spins b and c one obtains 

ip(b k+1 \X,a,d,Y\ +a 2 , -a B , +a 5 , -a 3 ) - tp{b k \ X, a, d,Y | +a 2 , -a 6 , +a 5 , -a 3 ) = , 

(70) 

Tp(c k+1 \X,a,d,Y\ +a 3 , -a 5 , +a 6 , -"2) - ip(c k \X,a,d,Y\ +a 3 , -a 5 , +a 6 , -a 2 ) = , 
where k = 1, 2, . . . , n — 1. 

4.3 Discrete non-linear equations 

In the limit ([4"T]) the partition function (f2"9"]) develops a typical quasi-classical asymptotics. Substi- 
tuting ([65]) into ([29]) and calculating the integral by the saddle point method one obtains, 

\ogZ = -^A({y}) + O(h ) (71) 

where the action 

A{y})= ^(a,6,c,d) (72) 

white stars 

where the sum over all white-centred stars and a, b, c, d denote the corresponding outer spins 
arranged as in Fig. HI Here {y} denotes an equilibrium configuration of spins on the black sub- 
lattice, which is determined by the Euler-Lagrange equations 

" A({y}) = , (73) 



5y(f 

where r is an integer two-dimensional vector numerating the internal sites on the black sub-lattice 
(remind that the spins are varied subject to the constraints (|10|) ) . Using the first expression from 
(|66|) one could bring these equations to the form 

— C ( uUy\a,b,c,d) =0 (74) 
by 

for every black site. The function C u ' v (y \ a, b, c, d) is defined in (|61|) and the symbols y, a, b, c, d 
denote the spins on a (black-centred) star, arranged as in Fig.[4]on the right side. The equations (I74j) 
are explicitly presented in ()62|) in an expanded form. Further, the definition ()66p for the Lagrangian 
£u,v(P'i b, c, d), entering the action ([72]) . involves the variational equations ([55]) . 

— £ ( u}v (x I a, b, c, d) =0 (75) 
ox 

for every white site. These equations are explicitly presented in (|59p . Combining the equations (|74p 
and ([75]) (respectively, their expanded forms ([59]) . ([62]) ). one gets a system of non-linear difference 
equations for all internal sites of the lattice (the boundary spins are fixed). Note, that each of these 
equations involves variables on five sites, belonging to a four-edge star. These classical integrable 
equations can be regarded as a generalization of the Laplace-type system associated with the Adler- 
Bobenko-Suris Q4 equation [13. 26j to the case of multi-component field variables. Their integrability 
is inherited from the integrability of the quantum model. It could also be independently established 
from the classical star-star relation (|63p (in particular, the latter implies the existence of an infinite 
set of local integrals of motion). 
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It appears that for Dirichlet (fixed) boundary conditions solutions of the classical equation 
([62]) possess the following property. We expect that in the bulk of a large lattice (i.e. away from 
the boundary) all solutions converge to some constant solution, independent of the boundary values 
of the spins. According to (|39p it is reasonable to expect that the action (|72p evaluated on this 
constant solution should vanish. The corresponding solution could be easily found. Assume that 
all spins on the lattice are equal to same vector 

x {c) = (xi,x 2 ,...,x n ), xj = - (j '•- —^—) > i = 1, 2, ...,n. (76) 



Then it is easy to check that all equations (|59j) . (|62|) are satisfied and that 

A a (x^ c \x^)=C{x^) = 0. (77) 

Therefore the action (j72|) exactly vanishes on this solution in complete agreement with ()39p . It 
would be interesting to verify whether this it is an absolute minimum of the action and, more 
generally, to study its global convexity properties. 

As an illustration consider the case n = 2. The two-component spin variables in this case 
contain only one independent continuous variable, so that it will be more convenient to use the 
scalar argument x, assuming that the corresponding two-component spin is equal to x = (—x,x), 
and similarly for all other spins. Eq. (|48p and f|51 1) simplify to 



*.(..») = -! P^log^^-ip^log**^), (78) 

#4(10 -\a\r) J n/2 <& A (w -ia\T) 



and 



7T 2 , 9 , ,7 77 



C( x ) = — + Ax 2 - 2ir\x\ , ~-<x<-. (79) 

Note that for n = 2 the function A a (x,y) becomes symmetric upon interchanging x and y. Using 
the above expressions in (I49D and (I66p one can bring the action (I72p to the form 

*£C(x s ) (80) 

(rs) s 

where the first sum is over all edges {rs), while the second sum is over all sites s. The variables a rs 
take one of the four values a.\, a 2 , ol$, 04, defined in (f56|) . depending on the type of the edge (rs) 
and the values of rapidity variables passing through the edge. The variational equations ()75p and 
for the white-centred stars, as in Fig. HI reduce to 



ip ai (x, a) + (p a2 (x, b) + ip a3 (x, c) + (p a4 (x, d) = 2n- 8x (81) 

where 

^ y) = 8l A„(x, ,) = i log M*-V + W*(* + » + «) (82) 

1 v^yx — y — \a.)v±\x + y — \a) 

Similarly, for the black-centred stars one gets 

¥>« 4 (y> °) + (y, &) + ^2 (V, c ) + ^ai (y> d) = 2vr - 8y (83) 
where the spins are arranged as in Fig. [H The constant solution (I76h in this case reads 

(84) 



-(c) _ 1 



4 
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for all sites of the lattice. Observing that 

A «G'i) =c (i) =o ' (85) 

one easily concludes that action (j80j) on this solution vanishes exactly A({xj}) = 0. Let us show 
that this is the minimum of the actiorJl. Expanding the latter in the vicinity of this solution 

x s = j + e s , e s ^0, (86) 

one gets 

AM) = (dflor.) " C 3 (ia rs )) {e r ~ e s ? + 2c ]T e 2 s + 0(e 3 ), (87) 

(rs) ' s 

where 

Ck(z) = T^-logtf^lr), fc = 3,4, (88) 
i oz 

and 

c = 2 + Ca(iai) + CsO^) + CaOas) + Cs(ia 4 ). (89) 
Remind that r is assumed to be purely imaginary, Imr > 0. Using the inequality 

2a 

<UQa) > > CsOa) > j-r , < a < n\r\/2 (90) 

TT\T\ 

and remembering that the sum of a's is constrained by (|57p . with % = — ivrr/2, one can check that 

c>0, 0<a fc <?7o, fc = 1,2,3,4, (91) 

and thus that (|87p defines a positive-definite quadratic form. The constant solution (|84p could 
be though as an analog of the isoradial solution of the Hirota equations describing planar circle 
patterns (TT| |3"6"]. 



5 Conclusion 

We formulated a new solvable edge-interaction model of statistical mechanics with multi-component 
continuous spin variables taking values on a circle. The Boltzmann weights are given by Eqs. (|12|) . 
(I13p . The weights satisfy the star-star relation (|26|) . which ensures the integrability of the model. 
Currently, we claim this relation as a conjecture. We have verified this relation in a few orders 
of perturbation theory in the temperature-like variables (see Appendix [B]) and made extensive 
numerical checks for n = 3,4. We expect that a complete proof could be obtained by a rather 
straightforward generalization of the results of |20H24j . devoted to discrete-spin models connected 
with the sl(n) algebra. 

The star-star relation (I26p implies the validity of the Yang-Baxter equations in the vertex (122p 
and IRF (j27|) forms. The partition function per edge of the lattice is calculated exactly in the 
thermodynamic limit by means of the inversion relations [33 35j . The result is included into the 
normalization of the Boltzmann weights given by (|14p , (|15p . With this normalization the bulk free 
energy of model vanishes, see (|39p . 

5 Since the expression (|79[1 is non-analytic for x = 0, one needs to check the case where all spins are vanishing 
a; s = 0. Careful considerations show that the value of the action for this configuration is greater than zero. 
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In Section S] we study the quasi-classical (or the low-temperature) limit when one of the elliptic 
nomes tends to unity, p — > 1, and obtain the classical version (Eqs. (|63p and (|68p ) of the star-star 
relation. 

Here we have only formulated our model on a simple square lattice, but the latter can be replaced 
with a rather general planar graph. The required construction is similar to that described in [13] with 
some modifications accommodating the two types of rapidity lines shown in Fig. [2j The partition 
function then will be invariant under certain deformations of the planar graph connected with the 
star-star (|26|) and inversion (I17D relations. This property is known as Baxter's Z-invariance |37j . 
Evidently, this property continues to hold in quasi-classical limit. As a result the action of the 
classical system (|72|) evaluated on the stationary configuration will be invariant under the "star- 
star" transformations connected with the classical star-star relation (|68p . These moves are similar 
to the "star-triangular" transformations for the classical action introduced in |10j . Remind that, 
geometrically, the classical star-triangle relation can be associated with a tetrahedron formed by 
a set of six face diagonals of a cube. By this reason the star-triangular moves are often referred 
to as "cubic flips" of faces of a quadrilateral surface [38] • From this point of view the star-star 
relation (I68j) could be associated with "rhombic-dodecahedral flips" of quadrilateral surfaces, since 
the octahedron in Fig. [5j representing the various terms in (|68p . forms a set of face diagonals of a 
rhombic dodecahedron. 

As mentioned in the introduction, for n = 2 the model provides a master solution [13] of the 
Yang-Baxter equation, which contains all known edge-interaction models with single-component 
spins as particular cases. It turns out that for n > 2 the model also possesses a similar "master" 
property. The Boltzmann weights (|12p allow a large number of interesting limiting cases. For 
example, one could consider a low-temperature limit when one of the nomes tends to a root of 
unity, 

p^e m/N , AT = 2,3,4,... . (92) 

This limit is similar, but more complicated than the quasi-classical limit considered in Section HJ 
In this case one obtains a "hybrid" model which couples a classical integrable system, involving 
continuous spin variables, and an Ising-type model of statistical mechanics with discrete multi- 
component spins variable taking values in (Zjy) ,l_1 - The Boltzmann weights are expressed with 
elliptic theta functions. In general, the emerging model is spatially inhomogeneous, since its Boltz- 
mann weights depend on solutions of the classical equation of motion H. In a particular homoge- 
neous case connected to a constant classical solution (similar to (I76f> ) this model reduces to the 
s/ n -generalization |24] of the Kashiwara-Miwa model [2j. In the trigonometric limit this elliptic 
model further reduces to the sZ n -generalized chiral Potts model |2ip22j. which is in its turn equiva- 
lent to the A r -state Zamolodchikov-Bazhanov-Baxter model [20, 39] on a cubic lattice (the number n 
becomes the size of the lattice in the "hidden" third dimension) . We hope to consider all these limits 
and connections in the future (the n = 2 case has been thoroughly studied in [H]). It would also 
be interesting to explore possible connections of our results to the discrete-spin models considered 
in 00]. 

6 These classical equation of motion are very similar to Eas. ()75p and (|74p . Conceptually the hybrid model described 
here is very much similar to a model of quantum field theory on a non-trivial classical background. 
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Appendix A. Inversion relation 

Theorem 11 from [16] (related to ^4 n _i root system) can be re- written in our notations as 

n 

dyW a (x,y)Wp(y,z) JJ $(0- - yj + ma) - a + in/3) = 
Kn(v — a — f3) $(ir/ — ma) &(\r] — in/3) 



X,Z 



$(ir/ - m(a + f}))K n (r] - a) K n (r] - /3) 



Appendix B. Series expansions of in powers of p and q 

Define power sums of exponents of the spin variables, entering the star weight ([2 



« = + 5(7-0), v = -5(7 + 0), 

«' = 5 + 0-5(7-0), ^ = +5(7 + 0) 



(A.l) 



x $(cr - a?j + ia) <&(zj - cr + i0) &(xj - a + i(n/3 - a)) <3?(cr - Zj + i(na - 0)) . 

Taking the limit /3 — >• — a and using (|16p and 

<&(*)<!>(-*) = 1 , Um $( 7 ~ ' n f = 1 , W (x,z) = - y2S(x,a(z)) , (A.2) 

one obtains the inversion relation (T[ 



^ fc = ^e 2i ^, B A = JV ifc6 ', * fc = ^e 2ifc % k = l,2,.... (B.l) 

Consider the physical regime ([2]), fl3Qf) . Then rapidity variables can be parametrized as 



(B.2) 



where the parameters s a, /3 and 7 are real. Moreover, assume 

I Rea| < 77, p ~ q ->■ 0, (B.3) 
Using the product expression ([J]) for the elliptic gamma-function it is not difficult to obtain the 
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following expansions for the Boltzmann weights (|12|) and (|13|) . 

« n ( J7 /2 + a)W„ /2+a (a,6) = 1 + \ (pq)* e 2a A\ B 1 + pq e 4a {(^) 2 (BJ 2 + ^ B 2 ] 

+ i (pq)l e 6a {(^) 3 (SJ 3 + 2^B 3 + S X S 2 } 
+ (pq)5 (p 2 + q 2 ) e 2Q 4? B 1 - (pq)i e" 2a ^ B{ 

+ ^ p 2 q 2 e 8a {(^) 4 (^) 4 + 6 (A{) 2 A\ (B x fB 2 + 3 (A* 2 ) 2 (B 2 ) 2 + 8AtA* 3 B,B 3 +6A* 4 B 4 ] 

+ pq (p 2 + q 2 ) e 4a (A*,) 2 (B,) 2 - pW At A, B\ B x + 0(p 5 ) , 

where the order of the correction terms here and below is shown for p ~ q, 

K n (a)W a (a,b) = l + pq(l + p 2 + q 2 ){e 2a ^S 1 -e~ 2a A 1 B$} - p 2 q 2 At A 1 B* 1 B 1 

+ p 2 q 2 { cosh 4a (At) 2 (B^ 2 + sinh4a A 2 B*\ + 0(p 6 ) , 



(B.4) 



(B.5) 



and 



7r n_1 n! §(x) H (2sm( Xj - x k )) 2 = 1 - (p 2 + q 2 ) {x?X 1 - l} + p 2 q 2 {x£X 1 - l} 2 

+ \ (p 4 + q 4 ) {4 -±X*X X + (X 1 *) 2 (X 1 ) 2 - X^X 2 } + 0(p 6 



j<k 



(B.6) 



(B.7) 



From (I15p one obtains 

k( V /2 + a) = 1 + (pq)? e 2nQ + 0(p n + 2 ) 

The definition of the IRF weight f)28|) requires an evaluation of integrals over the central spin in 
the star. The results depend on the value of n. Below we restrict ourselves to the case n = 3. For 
any function f(x) of the spin variable x define the following integral 

J[f(x)]= — x j [ ^JJ(2sin(xj -x k )) 2 ^j f(x)dx, dx = dx\ dx 2 ■ ■ ■ dx n -\ . (B.8) 

Below we will need the following integrals 

J[l] =j\x x X{ 



(B.9) 



J [(^(AT) 2 ] =J[X 2 X*}= 2, 
which are given for a general n and specific integrals 

J[(X x f] = J[(Xff\ = 1, 

J[X X X 2 \ = J[X*X 2 ] = -1, 
valid for n = 3. Using the above expansions one obtains for the IRF weight for n = 3 ([28 



(B.10) 



S(c)S(6) 



-I / a 6 



c 0? 



1 + pq P + (pq)2 Q 



(B.ll) 



p 2 q 2( 1+ l 7 ,2 + l^_ 5 A + l pq(p 2 +q 2 )p + 0(p5) 
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where we used the following notations 

V = (e 2i/3 Ai + e" 2 " 3 Pi) (e 2i ^ B\ + e" 2 ^ C{) + (e" 2i/3 A\ + e 2i/3 PJ) (e" 2 ^ B x + e 2i ^ Ci) 

Q = e 6Q (e- 2i/3 ^4iPJ + e 2i/3 ^4JPi) + e" 6 " (e^B^f + e-^BJCi) 

K = (e^ A 2 + e" 4i/3 V 2 ) (e 4i ^ 85 + e" 4 ^ C* 2 ) + (e" 4i/3 ^ + e 4i ^ D|) (e" 4 ^ £ 2 + e 4i ^ C 2 ) 

5 = (e 2i/3 A + e" 2i/3 Pi) (e" 2i/3 ^ + e 2i/3 PJ) + (e" 2 ^ B x + e 2i ^ C x ) (e 2i ^ £J + e" 2 ^ Cf ) 

(B.12) 

All additional integrals (beyond (|B.9|) and (|B.10p ). required for (|B.lip in the given order, exactly 
vanish. Note all terms in (IB. lip are manifestly real and thereby confirm the validity of the star-star 
relation (|26p in pertubation theory. 
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